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ABSTRACT 


With  the  increasing  use  of  bonded  repairs  to  restore  the  structural  integrity  of  ageing 
aircraft  the  question  arises  as  to  the  interaction  effects  when  repairs  are  located  close 
together.  Using  the  Finite  Element  (F.E.)  method  a  study  has  been  carried  out  for  the 
interaction  between  two  idealised  circular  repairs.  The  interaction  involves  the 
increase  of  the  sheet  stress  just  outside  the  patch.  It  has  been  found  that  the  tandem 
orientation,  with  respect  to  the  applied  load,  is  the  most  severe  configuration.  In  this 
case,  for  most  practical  repairs,  the  interaction  may  result  in  increases  of  the  sheet 
stress  by  40%  for  very  close  separation  distances.  It  has  also  been  found  that  certain 
combinations  of  bi-axial  load  can  also  significantly  influence  the  interaction  effect. 
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Analysis  of  the  Interaction  Effect  for  Bonded 

Repairs 


Executive  Summary 

The  repair  of  cracked  structures  using  bonded  repairs  has  resulted  in  considerable 
aircraft  lifetime  extension  and  hence  cost  savings.  The  long  fatigue  life  of  a  bonded 
repair  is  attributed  to  the  lack  of  rivets  and  hence  stress  concentrations  in  the  repair. 
However  the  widespread  use  of  these  repairs  for  ageing  aircraft  may  result  in  repairs 
being  located  close  together.  The  purpose  of  this  study  is  to  determine  the  factors  that 
influence  the  sheet  stress  just  outside  the  patch. 

The  finite  element  method  is  used  to  analyse  the  problem,  and  for  convenience  the 
reinforcement  has  been  modelled  to  be  a  circular  patch  with  constant  thickness.  The 
first  geometric  configuration  considered  is  the  tandem  configuration  with  respect  to 
the  applied  load.  This  is  the  most  severe  configuration  and  for  most  practical  repairs, 
the  interaction  may  result  in  increases  of  the  sheet  stress  by  40%  for  very  close 
separation  distances.  It  has  also  been  found  that  certain  combinations  of  bi-axial 
applied  load  can  also  significantly  influence  the  interaction  effect.  In  the  case  of  a 
periodic  array  of  repairs  in  the  tandem  configuration  under  a  uni-axial  load  the  highest 
sheet  stresses  are  increased  by  60%  for  close  separations.  In  the  case  of  the  side-by-side 
configuration  the  sheet  stresses  are  increased  only  by  19%  for  most  practical  repairs, 
and  correspond  to  large  separation  distances.  A  finite  width  effect  has  been 
investigated  for  the  tandem  configuration  and  it  has  been  found  that  for  small  widths 
the  sheet  stress  is  lowered  by  as  much  as  13%.  For  all  configurations  considered  over 
stiff  repairs  gave  the  maximum  stress  increases. 

The  design  of  bonded  repair  procedures  is  well  documented  in  the  RAAF  Engineering 
Standard,  C5033.  However  the  current  load  attraction  factor  used  for  bonded  repairs 
does  not  take  account  of  the  separation  distance  between  repairs.  The  purpose  of  the 
work  presented  here  is  to  present  design  curves  which,  when  used  together  with 
material  S-n  data,  will  enable  the  RAAF  to  more  accurately  assess  the  fatigue  life  of  the 
sheet  just  outside  repairs. 
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1.  Introduction 


The  repair  of  cracked  aircraft  structures  using  bonded  repairs  has  resulted  in 
considerable  aircraft  lifetime  extension  and  hence  cost  savings.  The  long  fatigue  life  of 
a  bonded  repair  is  attributed  to  the  lack  of  rivets  and  hence  stress  concentrations  in  the 
repair.  However  the  widespread  use  of  these  repairs  for  ageing  aircraft  may  result  in 
repairs  located  close  together.  The  purpose  of  this  study  is  to  quantify  the  interaction 
effect.  The  existence  of  this  effect  has  been  shown  in  Chow  and  Atluri  [1]  for  limited 
geometries.  In  the  work  carried  out  here,  a  parametric  study  is  carried  out.  In  this 
study  the  increase  in  the  stress  in  the  sheet  just  outside  the  patch  is  the  quantity  being 
considered.  The  orientation  of  the  repairs  is  considered  with  respect  to  the  applied 
load.  Furthermore  biaxial  loadings  will  be  considered.  The  finite  width  effect  is  also 
considered  for  two  different  configurations.  The  analysis  is  considered  as  linear  elastic. 
The  design  of  bonded  repair  procedures  is  well  documented  in  the  RAAF  Engineering 
Standard,  C5033  [2].  However  the  current  load  attraction  factor  used  for  bonded 
repairs  does  not  take  account  of  the  separation  distance  between  repairs.  The  purpose 
of  the  work  here  is  to  present  design  curves  which,  when  used  together  with  material 
S-n  data,  will  enable  the  RAAF  to  more  accurately  assess  the  fatigue  life  of  the  sheet 
just  outside  repairs. 


2.  F.E.  Methodology 

The  representation  of  the  repair  in  this  analysis  is  simplified  and  does  not  include  the 
adhesive,  in  accordance  with  stage  1  of  the  two-stage  design  analysis  proposed  by 
Rose  [3]  for  repair  design.  For  convenience  the  reinforcement  has  been  modelled  to  be 
a  circular  patch  with  constant  thickness  using  a  2D  linear  elastic  finite  element 
analysis.  Only  in-plane  loads  are  considered,  and  out-of-plane  secondary  bending  has 
been  restrained.  Thus  the  analysis  is  most  appropriate  for  two-sided  repairs,  or  for 
cases  where  out-of-plane  deflection  is  restrained  by  stiff  sub-structure  Rose  et.  al.  [4]; 
one-sided  repairs  introduce  geometrical  non-linearity  and  new  length-scales  which  are 
not  considered  here  Wang  et.  al.  [5].  Also  it  was  assumed  that  the  hole  or  crack  being 
reinforced  does  not  affect  the  load  being  drawn  into  the  patch.  As  a  result  a  simple 
model  has  been  developed  in  which  the  variables  are  firstly  the  ratio  of  the  combined 
Young's  modulus  of  the  patch  and  sheet  beneath,  E0,  to  the  sheet  modulus,  Es,  as 
indicated  in  Fig.  1,  and  secondly  by  the  separation  distance,  2 S,  between  the  two 
repairs  (Fig.  2).  From  equilibrium  considerations  the  E0  /  Es  ratio  for  any  repair,  as 
shown  in  Fig.  1,  is  given  by: 

Eo.  _  Ests+EPtp 

Es  Ests 
where 

Ep  Young's  modulus  for  the  patch 
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ts  thickness  of  the  sheet 

tp  thickness  of  the  patch 

For  most  repairs  Eptp  =  Ests ,  hence  E0  /  Es  =  2.  For  an  over-stiff  repair  Eptp  >  Ests. 


Figure  1.  The  geometry  for  an  equivalent  stiffness  representation. 

Load  cases  considered  are  for  remote  ex x  and  c  Y  stresses  and  a  combination  of  both 
a  x  and  aY  .  The  axis  system  used  is  shown  in  Fig.  2. 


▲ 


Figure  2.  Two  repairs  in  an  infinite  sheet:  axis  system  and  locations  of  interest. 
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The  separation  between  the  repairs  is  expressed  as  a  ratio  of  S/R.  Shown  in  Fig.  2.  is  a 
close  up  view  of  the  F.E.  mesh  for  a  separation  ratio  of  1.1.  A  range  of  separation  ratios 
of  1.1  to  5  will  be  considered. 


Figure  3.  Close  up  view  of  mesh  for  S/R=l.l 


3.  Results 


This  configuration  is  shown  in  Fig.  4  and  corresponds  to  an  applied  stress  of  crx  =  1 
and  <yy  =  0 .  The  stress  contours  indicate  that  the  maximum  principal  stress,  cTj , 
occurs  at  location  C  (see  Fig.  2).  Also  it  happens  that  the  maximum  principal  stress 
occurs  in  the  X  direction.  The  results  for  this  configuration  with  respect  to  the  applied 
load  are  shown  in  Fig.  5  and  indicate  a  strong  interaction  for  close  separation  and  high 
E0  /  Es  values,  and  maximum  stress  ratio  of  1.60  was  achieved.  Normally  a  repair 
would  correspond  to  E0  /  Es  =  2  and  hence  would  correspond  to  a  maximum  stress 
ratio  of  1.4.  At  large  separation  distances  cr,  approaches  an  asymptotic  value  which 
corresponds  to  that  of  a  single  repair. 

3.1  Tandem  configuration 


This  configuration  is  shown  in  Fig.  4  and  corresponds  to  an  applied  stress  of  a x  =  1 
and  (7 y  =  0 .  The  stress  contours  indicate  that  the  maximum  principal  stress,  cr, , 
occurs  at  location  C  (see  Fig.  2),  and  this  location  is  independent  of  the  S/R  ratio.  Also 
it  happens  that  the  maximum  principal  stress  occurs  in  the  X  direction.  The  results  for 
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this  configuration  with  respect  to  the  applied  load  are  shown  in  Fig.  5  and  indicate  a 
strong  interaction  for  close  separation  and  high  E0  /  Es  values,  and  maximum  stress 
ratio  of  1.60  was  achieved.  Normally  a  repair  would  correspond  to  E0  /  Es  =  2  and 
hence  would  correspond  to  a  maximum  stress  ratio  of  1.4.  At  large  separation 
distances  cr,  approaches  an  asymptotic  value  which  corresponds  to  that  of  a  single 
repair. 


Figure  4.  Tandem  configuration:  stress  interaction  betiveen  two  circular  repairs,  where  cr=crx-l 
and  S/R=1.2.  (A=.878,  B=.928,  C=.977,  D=1.027,  E=1.076,  F=1.126 ,  G=1.175,  H=2.225, 
1=1.274,  J=1.324) 


1.60 
1.55 
1.50 
1.45 
1.40 
b*  1.35 
tT  1.30 
1.25 
1.20 
1.15 
1.10 


Figure  5.  Tandem  configuration:  variation  of  stress  ratio  a\/ax  with  S/R  ratio. 
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or 

t 


I 

or 


Figure  6.  Side  by  side  configuration:  stress  interaction  behveen  two  circular  repairs,  ivhere 
<jx=l.  and  cTy=0.  and  S/R=1.2  (A=.810,  B=.848,  C=886,  D=.925,  E=.963,  F=1.001,  G=1.039, 
H=l. 078, 1=1.116,  J=1.154) 

3.2  Side  by  side  configuration 

This  configuration  is  shown  in  Fig.  6  and  corresponds  to  the  applied  loading  crx  =  0 
and  gy  -  1.  In  this  case  the  stress  contours  for  the  maximum  principal  stress  show 
that  the  maximum  value  occurs  at  location  B  (see  Fig.  2),  and  this  location  is 
independent  of  the  S/R  ratio.  Also  it  happens  that  the  maximum  principal  stress 
occurs  in  the  Y  direction.  The  results  for  a  study  of  the  interaction  are  shown  in  Fig.  7. 
These  results  are  different  to  the  previous  case  in  that  the  maximum  stress  ratios  occur 
at  large  separations,  due  to  shielding,  and  show  an  asymptotic  behaviour.  As  before, 
the  maximum  stress  concentration  in  the  sheet  occurs  for  high  E0  /  Es  values.  The 
maximum  stress  ratio  achieved  was  1.28.  Again,  for  most  repairs  (E0  /  Es  =  2.)  and 
the  maximum  stress  ratio  is  1.19. 
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Figure  7.  Side  by  side  configuration:  variation  of  stress  ratio  cr,  /  aY  with  S/R  ratio. 

3.3  Periodic  array 

This  is  the  representation  of  an  infinite  number  of  repairs  in  the  X  direction,  as  shown 
in  Fig.  8,  and  is  equivalent  to  the  tandem  configuration.  The  height  of  the  sheet  is  very 
large  in  comparison  to  the  radius,  and  hence  can  be  assumed  to  be  infinite.  The  loading 
is  carried  out  by  applying  uniform  displacements  in  the  X  direction  along  the  face  of 
the  sheet.  This  results  in  a  variation  of  a x  stress  in  the  Y  direction  and  from  this  an 
average  stress  has  been  calculated.  Results  for  a  range  of  Ea  /  Es  are  presented  in  Fig. 
9  and  correspond  to  point  A  shown  in  Fig.  8.  Although  the  separation  distance  is 
considered  in  terms  of  W/R  ratio's,  it  is  apparent  that  at  close  separations  maximum 
stress  ratio's  of  1.9  are  achieved  at  point  A,  with  a  value  of  E0  /  Es  =  3.0.  However  for 
most  repairs  the  maximum  stress  ratio  would  be  1.55. 
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Figure  8.  Periodic  array  in  x  direction  (tandem  configuration ),  infinite  in  y  direction. 


Figure  9.  Results  for  -periodic  array,  stress  ratio  oj/cw  for  W/R  ratios. 

3.4  Combination  of  loads 

A  bi-axial  stress  field  has  been  considered  in  which  the  applied  stresses  are  ctx  =  1  and 
ay  =  -1,  corresponding  to  the  geometry  in  Fig.  2.  These  results  are  shown  in  Fig.  10.  In 
this  case  the  maximum  interaction  occurs  at  location  C  (see  Fig.  2)  corresponding  to  the 
lower  S/R  ratios.  Again  the  repairs  with  large  Eo/Fs  ratios  result  in  the  largest  cr,  /  cr x 
ratios.  The  maximum  principal  stress  also  happens  to  occur  in  the  X  direction.  This 
combination  of  loads  results  in  a  stress  interaction  equal  to  that  of  the  periodic  array. 
For  the  case  where  a x  =  (Ty  =  1 ,  for  the  geometry  shown  in  Fig.  2,  results  obtained  are 
shown  in  Fig.  11.  In  comparison  to  Fig.  5  for  the  uni-axial  case,  the  combination  of 
stresses  alleviates  the  interaction  effect.  Clearly,  the  direction  of  the  load  components  is 
important. 
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Figure  10.  Bi-axial  loading:  S/R=l.l,  <JX  =  1,  aY  =  -1,  stress  ratio  at  point  C. 


Figure  11.  Bi-axial  loading ,  S/R=l.l,  ax  =  l,ay  stress  ratio  at  point  C. 

3.5  Finite  width  effect 

Consider  the  tandem  configuration  shown  in  Fig.  12.  In  this  case  the  S/R  ratio  has 
been  set  to  1.1  since  it  is  already  known  that  the  highest  interaction  occurs  with  this 
geometry.  The  results  shown  in  Fig.  13  show  a  percentage  change  of  sheet  stress 
(location  C)  with  a  variation  of  width,  W,  (expressed  as  a  ratio  W/R).  The  greatest 
reduction  in  sheet  stress  occurs  for  W/R  ratios  approaching  1,  and  for  higher  E0  /  Es 
ratios. 
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For  the  case  of  the  side  by  side  configuration  shown  in  Fig.  14  the  S/R  ratio  adopted 
was  3  since  the  highest  stress  occurs  at  large  separations.  The  results  for  the  percentage 
change  of  sheet  stress  (location  B)  versus  W/R  are  shown  in  Fig.  15.  Again  the  largest 
reductions  of  sheet  stress  occur  for  W/R  ratios  approaching  1,  and  for  higher  E0  /  Es 
ratios.  In  comparison  with  the  tandem  configuration  the  finite  width  effect  is  not  as 
strong.  The  results  may  be  surprising  that  a  smaller  width  should  decrease  the  sheet 
stress,  but  an  explanation  is  that  with  a  narrow  width  the  load  attraction  ability  is 
limited. 


Figure  12.  Finite  width  effect :  tandem  configuration. 


Figure  13.  Finite  xoidth  effect:  %  Change  in  sheet  stress  for  tandem  configuration. 
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<?Y 
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Figure  14.  Finite  width  effect:  side  by  side  configuration. 
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4.  Validation 


The  purpose  of  the  validation  is  to  verify  the  accuracy  to  which  stresses  can  be 
predicted  just  outside  a  repair  for  a  given  mesh  density.  In  particular,  the  mesh 
density  used  in  the  interaction  study.  Consider  a  circular  repair  on  a  circular  2D  sheet 
as  shown  in  Fig.  16.  The  sheet  is  subject  to  a  uniform  pressure  load  applied  to  the 
perimeter  of  the  sheet  in  a  radial  direction.  The  patch  region  has  a  radius  r=Ri ,  while 
the  radius  of  the  circular  sheet  is  r=Ro  ■  Quantities  Ri ,  Eo  Jo  Mo  and  c 0  are  relevant  to 
the  patch  region  for  r  <  R0 , while  quantities  Ro  ,  Es  Js  Ms  and  as  correspond  to  the 
sheet  region  where  R;  <r  <R0. 


Plate  subject  to  radial  in-plane  stress  <j(r  =  R0) 


E  s  i 


E  oJ 


o  ,  o 


(b)  Side  view 

Figure  16.  Variables  for  closed  form  solution  for  a  circidar  repair  on  a  circular  sheet. 
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Firstly  the  radial  stresses  and  displacements  will  be  obtained  for  the  sheet  and  patch. 
Consider  the  configuration  with  properties  as  shown  in  Fig.  16.  The  governing 
equation  is  given  by  Timoshenko  &  Goodier  [6]  as: 


d2u  1  du  u 
dr 2  r  dr  r2 


(2) 


where  u  is  the  radial  displacement 
r  is  the  radial  distance 


The  solution  of  this  equation  is  given  by: 
for  r  >  R, 

C, 

u,  =  C2r  +  —  (3) 

r 


for  r  <  Rj 


uQ  =  Cxr 

where  C, ,  C2  and  C3  are  constants  of  integration. 

It  can  also  be  shown  from  [6]  that  the  radial  stress  is  given  by: 
for  r  >  R, 


<y  c  = 


(i-R  )L 


C3(l-  v) 
C2(  1  +  v)-  3~—  J- 


(4) 


(5) 


for  r  <  Rj 

^=^T^yc>(1+v) 

It  is  now  necessary  to  solve  for  the  three  constants  of  integration.  The  solution  of  these 
equations  must  satisfy  the  following  conditions: 

(a)  The  displacement  us  and  u0  is  equal  at  r  =  R;,  hence  from  equations  (3)  and  (4): 

a 

C2r  +  —  =  C^Rj  (7) 

r 

(b)  Equilibrium  must  be  maintained  across  the  boundary  at  r  =  R,  (note  that  ts  =  tQ ), 
using  equations  (5)  and  (6)  leads  to: 
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Es^s 


(1-v2) 


C2  (1  +  v) 


C3(l-v) 

R) 


E0t0 

(1-v) 


Q 


(8) 


(c)  Also  we  have  the  from  equ(6)  that  at  r  <  Rj : 

<7»  =  (ir^c'(1+,/)  (9) 

qqg  gta.ge  we  have  enough  information  for  the  evaluation  of  the  constants  C) ,  C2 
and  C3 ,  hence: 


(1-v) 

E0 


Q  = 


(l  -  v)cr0 


1  (1+!/) 

iY 

E0  2 

<E0 

~eJ 

Using  equs(5,  6)  for  cr,  and  <J0  we  obtain  the  following: 


(10) 


(11) 


(12) 


E* 

\  —  +  0.5 

f  Es) 
1-TT 

- 1 

+ 

.+ 

pO 

'(l-v) 

Eq 

V  Eq) 

L 

l  r  J 

(13) 


Results  for  two  Es  /  E0  ratios  are  shown  in  Fig.  18  together  with  the  F.E.  results.  The 
F.E.  mesh  is  shown  in  Fig.  17.  Very  good  agreement  between  the  closed  form  solution 
and  F.E.  results  has  been  achieved.  It  is  evident  that  the  stress,  just  outside  the  patch, 
can  be  predicted  to  a  high  degree  of  accuracy.  Since  a  similiar  radial  mesh  density  has 
been  used  in  the  interaction  study,  this  gives  confidence  in  the  results  that  have  been 
obtained  in  the  repair  interaction. 
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Figure  17.  F.E.  mesh  for  validation  of  circular  repair  on  a  circular  plate ,  Ri=162.mm, 
Ro=500.mm. 
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Figure  18.  Validation  of  F.E.  work  using  closed  form  solution  for  a  single  repair. 
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5.  Conclusion 


Firstly,  the  validation  using  a  closed  form  solution  for  a  single  repair,  also  derived 
here,  gives  very  good  agreement  with  the  F.E.  results.  This  gives  confidence  that  the 
stress  predictions  in  the  sheet,  just  outside  interacting  repairs,  is  accurate.  Results  for 
the  repair  interaction  indicate  that  maximum  stress  increases,  depend  on  both  the 
orientation  of  the  repair  with  respect  to  the  applied  load,  and  separation  distance 
between  repairs.  The  side  by  side  configuration  with  respect  to  the  applied  load  gave 
stress  increases  of  between  12-27%,  where  the  maximum  values  occurred  at  the  larger 
separation  distances.  In  the  case  of  the  tandem  configuration  with  respect  to  the 
applied  load,  stress  increases  of  between  23-60%  were  obtained  and  the  maximum 
correspond  to  small  separation  distances.  In  both  orientations  over  stiff  repairs  gave 
the  maximum  stress  increases.  Thus,  shielding  and  load-shedding  behaviour  for 
neighbouring  repairs  is  in  complete  contrast  to  that  for  neighbouring  holes  or  cutouts, 
which  could  be  regarded  as  limiting  cases  for  which  Eo/Es— >0,  whereas  for  repairs 
Eo/Es  >  1. 

Bi-axial  loading  combinations  can  either  increase  or  decrease  the  stress  field 
interactions.  If  the  applied  component  stresses  are  of  the  same  sign  interaction  stresses 
are  decreased;  while  applied  stresses  of  opposite  signs  result  in  increased  interaction 
stresses.  Results  from  the  study  of  the  periodic  array  in  one  direction  have  shown  that 
this  configuration  gives  the  highest  stress  increases  for  uni-axial  loading,  and  also  is  a 
tandem  orientation  with  respect  to  the  load.  Furthermore  it  has  been  shown  that  the 
location  of  the  maximum  stresses  are  independent  of  the  S/R  ratio.  An  investigation  of 
the  finite  width  effect  has  shown  that  repairs  in  sheets  with  small  widths  result  in 
smaller  stresses  than  in  large  width  sheets.  In  this  case  the  repair  is  no  longer  able  to 
attract  additional  load  from  the  surrounding  structure. 

The  design  of  bonded  repair  procedures  is  well  documented  in  the  RAAF  Engineering 
Standard  C5033.  However  the  current  load  attraction  factor  used  is  independent  of  the 
separation  distance  between  repairs.  The  design  curves  presented  here,  when  used 
together  with  material  S-n  data,  will  enable  the  RAAF  to  more  accurately  assess  the 
fatigue  life  of  the  sheet  just  outside  repairs. 
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